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Abstract 

The Dirac structure of confinement is shown to be of timelike-vector nature 
in the heavy quark limit of QCD. This stands in contradiction with the phe- 
nomenological success of the Dirac scalar confining potential. A resolution 
is achieved through the demonstration that an effective scalar interaction is 
dynamically generated by nonperturbative mixing between ordinary and hy- 
brid QQ states. The resolution depends crucially on the collective nature of 
the gluonic degrees of freedom. This implies that dynamical gluonic effects 
are vital when attempting to incorporate fine structure in models of the QQ 
interaction. 
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I. INTRODUCTION 



Although it has been postulated for more than 30 years, the phenomenon of quark 
confinement remains an enigmatic feature of QCD. Quenched lattice gauge theory and heavy 
quark phenomenology indicate that the static {m q » Aqcd) long range potential should 
be linear with a slope of b « 0.18 GeV 2 . The order 1/m 2 , quark-antiquark long range spin- 
dependent (SD) structure has also been studied. Comparison with spin splittings in the J/ip 
and T spectra |TJ] indicate that the spin-dependence can only arise from the nonrelativistic 
reduction of a scalar current quark-antiquark interaction. This picture is also supported by 
calculations of the long range spin- dependent effective potentials on the lattice [§]. 

Unfortunately, little analytical progress has been made on this problem. The framework 
for most investigations on this subject was provided by the Wilson loop approach of Eichten 
and Feinberg [f| who extended the analysis of the spin-independent potentials by Brown 
and Weisberger []|]. The standard parameterization for the long range SD quark-antiquark 
interaction introduced in Ref. [Q is given by 

V SD {r) = ( (Tq ' 1Lq - a ~ q ' Lq ] (-— + -—] + ( °~ q ' Lq - ^"M ( ldV A 
\ Arrig 4m? J \r dr r dr J \ 2m q mq 2m q m q - ) \r dr ) 

+ i2^S 3<Tq ' iaq ' * " aq ' CT ^ 3(r) + T^" 9 ' "^ 4(r) - (1) 

Here e = e(r) is the static potential, r = |r| = |r g — r^| is the QQ separation and the 
Vi = Vi(r) are determined by electric and magnetic field insertions on quark lines in the 
Wilson loop expectation value (explicit expressions are given in Section II. C). As shown by 
Gromes covariance under Lorentz transformation leads to a constraint exists between the 
SD potentials, 

e(r) = V 2 (r)-V 1 (r). (2) 

Other more or less fundamental relations were also derived IggQ. 

In a model approach, QQ interactions are typically derived from a nonrelativistic re- 
duction of a relativistic current-current interaction. As far as long range potentials are 



concerned only time-like vector or scalar currents are relevant ||. Performing a nonrela- 
tivistic reduction of a vector-vector interaction yields V\ — 0, V2 — e, V3 = e'/r — e", and 
V4 = 2V 2 e. Alternatively, the reduction of a scalar interaction yields 

V 1 = -e, V 2 = V 3 = V 4 = 0. (3) 

It is the alternation in sign of the combination V\ + V2 between vector and scalar currents 
which, through the analysis of the heavy quarkonia spectrum, enabled Schnitzer to identify 
the scalar interaction as the likely structure for confinement ||. 

Calculations [|K],[ll]] of the SD potentials based on sophisticated models of the Wilson 
loop typically yield results which are in agreement with Eq. (|3]). In particular the "minimal 
area law" model for the Wilson loop, a simple extension of the strong coupling limit of lattice 
QCD, leads to a picture which is very close to the classical flux tube model of Buchmuller 
fT2|| . He noted that if one assumes that the chromoelectric field is confined within a tube 
between the QQ pair then the magnetic field generated by the flux tube vanishes in the 
individual quark rest frames and the only contribution to the fine structure comes from the 
kinematical Thomas precession. This results in a SD structure which is identical to that of 
Eq. ©. 

A consistent picture of a QCD-generated effective scalar confinement interaction appears 
to be emerging. It is therefore disconcerting that attempts to build Hamiltonian-based 
models of QCD |13| seem to require vector confinement. For example, we have found |14| . |Tr) 



that it is impossible to construct a stable BCS-like vacuum of QCD when scalar confinement 
is assumed. This is problematical if one wishes to dynamically generate constituent quark 
masses. Furthermore, attempts at modeling chiral pions will be hindered by the explicit 
lack of chiral symmetry in a scalar interaction. These observations appear to stand in 
contradiction to the well established scalar confinement hypothesis. 

In the following we resolve this issue by first performing a Foldy-Wouthuysen reduction of 
the full Coulomb gauge Hamiltonian of QCD. This immediately establishes that the Dirac 
structure of confinement for heavy quarks is of a timelike-vector nature. It is commonly 



stated in the literature that vector confinement is ruled out by its spin-dependent structure. 
We wish to stress that one must be careful in this judgement. In particular, the spin- 
dependent interactions which are generated by QCD are more complicated than those given 
by the simple nonrelativistic reduction of an effective long range interaction. Indeed, we shall 
demonstrate that the scalar character of the spin splittings in heavy quarkonia is dynamically 
generated through effective interactions which crucially depend on the collective nature of 
the gluonic degrees of freedom. 

II. HEAVY QUARK EXPANSION OF H QCD 
Our starting point is the Coulomb gauge QCD Hamiltonian 

Hqcd = J dx^(x) [-ia ■ V + (3m] V>(x) + H YM + V c -g J dxtf(pc)a ■ A(x)^(x) (4) 
where 

V c = \g 2 f d*dyJ~ l p a (K)V a \^y;A)Jp\y), (5) 
K6(x, y; A) = (x, a|(V • P)- 1 (-V 2 )(V • V)~ l \y, b), (6) 



and 



l - J dx [j-'U^JU^) + B 2 (x)] . (7) 



Hym 

The degrees of freedom are the transverse gluon field A = A a T a , its conjugate momentum 
II = II a T a , and the quark field in the Coulomb gauge. The Faddeev- Popov determinant is 
written as J = Det[V • V], with V ab = \75 ab — gf abc A c being the covariant derivative in the 
adjoint representation, and the magnetic field is given by B a = V x A a + gf abc A b x A c . The 
static interaction Vc is the nonabelian analog of the Coulomb potential. It involves the full 
QCD color charge density which has both quark and gluon components, 

p a (x) = V t (x)T>(x) + f abc A b {x) ■ IT(x). (8) 



The most salient feature of the Coulomb gauge Hamiltonian is that all of the degrees of 
freedom are physical. This makes it especially useful for identifying the physical mechanisms 
which drive the spin splittings in heavy quarkonia. 



A. The Foldy-Wouthuysen Transformation 

We proceed by performing a Foldy-Wouthuysen transformation on the QCD Hamilto- 
nian. This is done in complete analogy to the quantum mechanical case where an operator 
is constructed which removes the interactions between upper and lower components of the 



quark wave function order by order in the inverse quark mass fL? |, except that the unitary 



transformation is now constructed in Fock space. The resulting Hamiltonian is given by 

Hqcd -> H FW = Jdx (m g /i t (x)/i(x) - m 5 x t (x)x(x)) + H YM + V c + H x + H 2 + . . . , (9) 
#i=2^/ rfx/lt ( x ) (D 2 - 9<r ■ B) h(x) -{h^ X ;m q ^ m q ), (10) 
H 2 = -^— f dxh?{x.)g<r ■ [E, xD]/i(x) + (h -> x, ™g -> rn q ). (11) 

In this expression h = (l+/?)^>/2 and x — (1 — f3)4>/2 denote the upper and lower components 
of the quark wave function and correspond to the annihilation and creation operators of the 
heavy quark and antiquark respectively. The ellipsis denotes terms which are either of 
0(l/m 3 ) or are spin-independent at order 1/m 2 . Finally D = zV + gA is the covariant 
derivative in the fundamental representation. The electric field contains both transverse 
and longitudinal components, E a = — n a + Ejj, where 

Ejj = -VA% - gVV~ 2 f abc A b ■ VA C (12) 

and 

A a Q ^)=g J dyV ab {^y-A)p\y). (13) 
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B. The Static Potential 



To leading order in the quark mass the Hamiltonian describes two static, noninteracting 
quarks. At O(m ) the Hamiltonian reduces to H = Hym + Vc- The eigenstates of H 
may be labeled by the quark and antiquark coordinates and by an index which classifies the 
adiabatic state of the gluonic degrees of freedom, n r , 

H \n r ; r q v q ) = e n (r)\n r ; r g r q ). (14) 

Note that we have made explicit the dependence of the gluonic degrees of freedom on the 
position of the quarks, r. 

The corresponding eigenenergies, e n (r) may be identified with the Wilson loop potentials 
calculated on the lattice. Thus for example, €o(r) is the Coulomb plus linear potential seen 



long ago [0|. Static hybrid states are collectively denoted \n r ; r g rg) with n r 0. In recent 



studies fl9| the lowest lying adiabatic hybrid potential, ex(r) has been evaluated. 

While both Hym and Vc may contribute to the linearly rising potential energy seen 
on the lattice, it is clear that the quarks may only interact with the flux tube via the 
nonabelian Coulomb interaction. Therefore the Dirac structure of confinement corresponds 
to 70 <8> 70 from the product of color charge densities (see Eqs. @ and (§)). As stressed in 
the Introduction, this appears to be at odds with 20 years of quark model phenomenology. 
Since the phenomenology is based on spin splittings, it will be instructive to examine the 
1/m 2 perturbative corrections to the static potential. 

C. Spin-dependent Potentials 

The spin-dependent first order correction to the static potential is given by 

fcn\ r ) = ^-^( n r; r q r q\ / d*h j (x)(T ■ [E, xD]/l(x)|?l n r q Tg) + (ft -> X', ™<q ~> m q) ■ (15) 

The order 1/m term is not considered because D 2 is not spin-dependent and the matrix 
element of cr • B vanishes. Eq. ([15]) may be simplified considerably as follows. Since we are 
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interested in spin- dependent terms only, the covariant derivative may be replaced with the 
ordinary derivative, zV, and the electric field may be replaced by — VA°. Contracting the 
fermion field operators and using Eq. (O) yields 



fc? } (r) ~ -^e^ k al(n r \Vi a (T a V a \r q , r q ; A)T b )\n r )V k ra + (/»-*; m, -> m q ; a q -> -a, 
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(16) 



where the matrix elements are over gluonic degrees of freedom only. The approximation sign 
is meant to serve as a reminder that the equality holds for spin dependent terms only. The 
above equation may be simplified further by using the following relation 

(n r \(Vi q g 2 T a V ab (r q ,r q ; A)T b )\n r ) = -V^e n (r). (17) 

The physical content of this relationship is simply the statement that Hym does not explicitly 
depend on the quark positions. The minus sign is due to contracting the antiquark operators. 
Quark line contraction also yields tadpole terms which vanish in the color singlet background 
and self energy diagrams which are subsumed into the leading order, 0(m), Hamiltonian. 
Making the appropriate substitution yields the standard classical plus Thomas precession 
spin orbit interaction 

= ( ^fnr - -^r- ( 18 ) 

\ 4m 2 , 4m? J r dr 

Thus the first term in Eq. (|l|), generalized to any adiabatic potential and therefore true for 

both ordinary and hybrid QQ states, is reproduced. At second order in perturbation theory 

the SD corrections are given by 

5e ( 2)(r) = IKir^-l^lm ;r r g -)| 2 ^ 
r^n en(r)-e m (r) 

In this expression there are two terms which correspond to the application of the magnetic 
field operator twice on a single quark or antiquark line. These matrix elements are not 
spin dependent since the product of the Pauli matrices collapses to unity plus a single Pauli 
matrix. Alternatively, the case where the magnetic fields act on different quark lines is 
nontrivial and yields 



a 2 

x,(2)i _ y i j v 
d6n lBB ~ 4m q m, a ^ X 

( y> (nrj r g r g | /rfxfo t (x).B , (x)fo(x)|m r ; r g r g ) (m r ; r g r g | Jrfyx^y^-^yMy)!^; r g r g ) 
Wn e n (r)-e m (r) 

(/i<-x)Y (20) 

If we define (g 2 times) the term in brackets as (f V J — \5 ij )V z + Id^VA then the third 
and fourth terms in the expression of the spin-dependent potential in Eq. (|1|) follow. If 
time-dependent fields are considered Eq. ( [ZUD may be rewritten as follows 

r T / 2 ,. ,., / 1 

-T/2 



f ' ( 1 \ 

lim g 2 / dtdt' — (n r ; r q r q \B % (r q , t)B J (r q , t')\n r ; r g r g ). 

^oo 7 -T/2 \ i / 

(21) 



Setting n r = in this expression yields a result which agrees with that of Eichten and 
FeinbergQ. 

There are four terms in Eq. (|19D which contribute to V\. These involve the application 

of the magnetic field and the covariant derivative on the same quark line, 

-9 ( y- (n r ; r g r g -| / cix/i t (x)D 2 /i(x)|m r ; r g r g ) (m r ; r g r g | / dyh)(y)(T ■ B(y)fo(y)|W r ; r g r g ) 
4m2^ n Cn ( r )_ em ( r ) 

#.c. + (/i -> x; ™ q -> m ? )J . (22) 
Using the following relationship (which holds for the spin dependent terms only) 

-^-/i t (x)D 2 /i(x) ~ -2^/i t (x)E(x) ■ V/i(x), (23) 



allows one to replace the first factor in Eq. ( p22|) with one involving the electric field. Per- 
forming the time integral and contracting the fermion field operators results in 
io\, Q 2 „• I v (n r \E l (r a )\m r ) (mJBi (r„)\n r ) \ , 

(24) 



1 Their expressions are in terms of Wilson loops and therefore project onto the ground state as T 
goes to infinity 
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If g 2 times the term in brackets is defined as 

dr r 

then this yields a spin orbit interaction in agreement with the second term in Eq. ([I]). 

Alternatively, if the time integral is retained after substituting Eq. ( |23| ) in Eq. (22) and 

the sum over energy denominators is also represented by another time integral, then the 

expression for V\ can be cast into the following form 

dVi ig 2 r T / 2 ft' — t\ 

*~dr = Sol J-T/2 dtdt ' \ ( n ^ r <7 r ?l B ( r 9>*) x E(r q ,t')\n r ;r q rq) (26) 

which agrees with the result of Eichten and Feinberg || . 

The remaining contributions to the second order energy correction correspond to an 
interaction between the quark and antiquark with both electric and magnetic operators at 
the vertices. Similar manipulations as for V\ yield the portion of Vsd proportional to V2 
once the following definition is made: 

which is equivalent to 

= lkn ^ f^dtdt' r-jr^j (n r ;r q r q \B(r q ,t) x E(r„ t')\n r ; r q r q ). (28) 

The Faddeev-Popov determinants have not been carried through the calculations shown 
above. They may easily be restored without changing any of the results. The time integral 
representations of the potentials shall prove convenient for subsequent calculations; however 
they are rather opaque. Alternatively, the application of the Foldy-Wouthuysen transforma- 
tion to the Coulomb gauge QCD Hamiltonian shows that these may be simply interpreted as 
nonperturbative mixing with hybrid states. This makes it clear that it is possible for non- 
perturbative dynamical physics to generate an effective spin-dependent interaction which 
mimics a scalar interaction. Actually demonstrating this requires that the matrix elements 
be evaluated. In the next section we propose to do this with the Flux Tube Model of Isgur 
and Paton. 
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III. MODEL EVALUATION OF THE SPIN-DEPENDENT POTENTIALS 



Before proceeding to a model evaluation of the matrix elements we note that it is possible 
to make some general statements on their expected properties. The results of lattice gauge 
theory make it clear that a flux tube-like configuration of glue exists between static quarks. 
If one thinks of this as a localized object with an infinite number of degrees of freedom, 
then it is apparent that V2 must evaluate to zero. This is because the electric field operator 
creates a local excitation in the flux tube at position r q (cf., Eq. ( ffip ). This must then 
be de-excited at r 5 by the magnetic field operator. However, the two operators become 
decorrelated because infinitely many degrees of freedom intervene. Similarly, the long range 
portions of V3 and V4 both vanish. Thus, by Gromes' relation (Eq. @), the only nonzero 
long range interaction must be given by V\ = — e. This is precisely the situation required for 
"scalar" confinement. It is therefore entirely plausible that an effective scalar confinement 
is generated by nonperturbative mixing with hybrids. Furthermore, the structure of the 
spin-dependent terms depends crucially on the nature of the ground state gluonic degrees 
of freedom and clearly favors a collective rather than a single particle picture of them. 

These simple expectations are borne out in an explicit model calculation. We shall 



employ the Flux Tube Model of Isgur and Paton [ 2DJ for this purpose. The model is extracted 
from the strong coupling limit of the QCD lattice Hamiltonian. The authors first split the 
Hamiltonian into blocks of distinct "topologies" (in reference to the possible gauge invariant 
flux tube configurations) and then make adiabatic and small oscillation approximations of 
the flux tube dynamics to arrive at an N-body discrete string-like model Hamiltonian for 
gluonic degrees of freedom. This is meant to be operative at intermediate scales a ~ 
where the strong coupling is order unity. The lattice spacing is denoted by a and there are 
N "beads" (or links) evenly spaced between the QQ pair. These considerations led Isgur 
and Paton to write the following model Hamiltonian 

N 1 h n N+1 

Hft = bo r+Y: ^-Pl + -f E (X(n) ~ x(n - l)) 2 (29) 

n=l 2b a 1 n=l 
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Here 60 is the bare string tension, p n is the conjugate momentum, and x( n ) is the transverse 
displacement vector at site n. It has been generalized to include color degrees of freedom: 
X — X c \( n )> c = 1 ■ ■ - 8, n = 1 ■ ■ ■ N, and A = 1,2. The ends of the string are fixed at n = 
and n = N + 1 by the static quark and antiquark positions and r = a(N + 1). 
In normal coordinates the Flux Tube Hamiltonian is 

N 2 ( 1 \ 
H FT = b r + ^J—^L + h ^uJ 2 n q 2 nX (30) 

where the normal mode frequencies are given by u n = | sin( 2( .^ 1 - ) ). This Hamiltonian may 
be trivially diagonalized with the following canonical transformation: 



a " A "V 2 g ™ A+ v^M^' (31) 



and one obtains 



1 7T 

Hft = u n al x a nX + br — + . . . (32) 

n\ a VlT 

The string tension has been renormalized by the zero point energy which also introduces a 
(divergent) constant and terms higher order in 1/r as shown above. 

Evaluating the spin dependent potentials in the Flux Tube Model requires explicit ex- 
pressions for the electric and magnetic fields in a flux tube. It is therefore necessary to extend 
the Flux Tube Model somewhat. We shall find it convenient to work in an "intermediate" 
coupling regime where we shall use the strong coupling lattice Hamiltonian for intuition 
and the weak coupling limit for the identification of the electric and magnetic fields. In the 
strong coupling limit the electric field operator simply counts links. It is therefore natural 
to map it onto a link displacement, 

E a x{n)^^{ X a x {n+l)-xl{n)). (33) 

For the sake of clarity, we shall take r q = and r 5 = rz from now on. This implies that 
A = x, y. The factor k is an arbitrary constant and will be identified later. 
The commutation relation between the electric and magnetic fields is 
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[E? (x), flj(y)] = ze^V^(x - y)(5 ab + (34) 
which implies that the magnetic field operator may be defined as the momentum conjugate 

to xx m 

B a Jn) = ——— — —. — -. (35) 

The e lJ of the commutator in Eq. fl34|) is taken into account by the index A, Xs{ n ) — ~Xy{ n ) 
and Xy{ n ) = Xx{ n )- Note that this relationship is physically sensible because the magnetic 
field operator maps onto the plaquette in the strong coupling limit and the application of a 
plaquette to a flux tube has the effect of moving a link one unit in the one of the directions 
transverse to the QQ axis with the magnetic field pointing in the other transverse direction. 
Substituting Eqs. (|33"D and (|35]) into Eq. (29) and setting k = ay/b~o yields 

H FT = b r+ l -jd* ((^(x)) 2 + (B a ± (x)) 2 ) . (36) 

This is reminiscent of the Coulomb gauge Hamiltonian where Hym involves transverse fields 
only and where the nonabelian Coulomb interaction has been replaced with a linear potential. 
Thus a satisfying consistency has been achieved with this approach. 

It is now a simple matter to write the field operators in terms of the collective phonon 
operators, a. One obtains 



(37) 



and 



Bt(xn, *) = -Vt £ sin(-^-n)v^ (< A e^ - o^e^-') (38) 

"I'm IV -T 1 



Substituting these expressions into Eqs. fl26|) and (Pq) , taking the phonon operator 
matrix elements in the 0-phonon number ground state, doing the time integrals, and 
evaluating the sum over modes yields 
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Vi(r) 



2a 2 



C F r 



(39) 



and 



V 2 (r) = lim AtC f — . 
2V ; Jv-oo2a 2 iV 



(40) 



In the strong coupling limit one has b = g 2 / (2a 2 )Cp so that the anticipated expression for 
V\ emerges in a natural way. Furthermore, Vi approaches zero like l/N; this is also true 
for V3 and V4. The latter point is illustrated in Fig. 1 where the correlation of electric 
and magnetic fields versus separation along the flux tube is shown. As expected the fields 
become completely decorrelated as the number of intervening degrees of freedom becomes 
large. Notice that this implies that a constituent gluon model of hybrids would not have 
been able to produce an effective scalar interaction. 

FIGURES 
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Fig. 1. Field Correlation Functions (for N=100). 

IV. CONCLUSIONS 

Spin splittings and lattice calculations indicate that confinement is scalar in nature. This 
conflicts with many relativistic models of QCD which require vector confinement. For exam- 
ple, a chirally symmetric interaction is needed if pseudo-Goldstone pions and spontaneous 
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chiral symmetry breaking are to be generated dynamically. Furthermore it appears to be 
impossible to build a stable vacuum with a scalar kernel. We have examined this issue with 
the heavy quark limit of the Coulomb gauge QCD Hamiltonian. This approach is physically 
intuitive and is simpler to interpret and implement than methods based on the Wilson loop. 
We found that the static confinement potential must indeed be a Dirac timelike vector. 
Effective scalar interactions are generated at order 1/m 2 by nonperturbative mixing with 
hybrid states. 

We have argued that the long range spin-spin (V3 and V4) and the vector-like spin-orbit 
potentials (V2) should all be zero since they involve field correlation functions evaluated 
between quark and antiquark. This statement follows by assuming that the gluonic degrees 
of freedom collapse into a flux tube-like configuration, as shown by lattice gauge theory. 
Alternatively, the scalar-like spin-orbit potential (Vi) is proportional to the matrix element 
of the electric and magnetic fields evaluated at the same point and hence is expected to 
be nonzero. Explicit calculations of the relevant matrix elements were carried out in the 
Flux Tube Model. The model was extended to include color degrees of freedom and to map 
the chromoelectric and chromomagnetic fields to flux tube phonon operators. The results 
obtained were in agreement with our general arguments and with Gromes' relation. 

A consistent picture of the Dirac structure of confinement has emerged. The static central 
potential is timelike vector while the spin-dependent structure mimics the nonrelativistic 
reduction of an effective scalar interaction. This implies that it is incorrect to employ a scalar 
confinement kernel when doing calculations with light quarks. Note however that it would be 
acceptable to use scalar confinement when working explicitly in the chiral symmetry broken 
phase, i.e., with constituent quarks in the nonrelativistic limit. The work presented here 
also implies that a constituent gluon picture of hybrids will yield incorrect results for certain 
observables. For example, V\ and V 2 would be of comparable magnitude in a constituent 
gluon model. In general, these types of models must fail when nonlocal properties of the 
gluonic configuration are considered. Alternatively, it is possible that they perform quite 
well when evaluating global properties of gluonics such as the hybrid spectrum. 
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The application of these results to spin splittings in heavy quarkonia is not straight- 
forward. For example, there is the possibility of large 0(1/ m 3 ) corrections to the split- 
tings. Light quark loop effects may also contribute to spin-dependent forces. It is, unfortu- 
nately, rather difficult to quantify these effects. Perhaps the best hope is with high precision 
NRQCD lattice calculations. Light quark loops may be studied by examining shifts in the 
predicted spin splittings between quenched and unquenched calculations. One naively ex- 
pects that these effects will be independent of the (heavy) valence quark mass. Alternatively, 
the higher order corrections should become larger as one moves from the T system to the J/ if) 
system (the quenched approximation may used to test this). Lattice results for the quenched 
T spectrum [f22| appear to be very close to experiment - supporting the idea that light quark 



loop effects are negligible (or at least may be absorbed into the parameters). Indeed, the 
calculated value of the ratio of P-wave splittings p = [M( 3 P 2 ) - M ( 3 Pi)]/[M( 3 Px) - M( 3 P )] 
yields 0.7(3) |23fl , in good agreement with the experimental value of 0.66(2). Unfortunately, 
this may not be regarded as definitive since p was evaluated at (3 = Q/g 2 = 6.0; a similar 
calculation at (3 = 5.7 gives p = 1.4(4), very far from the experimental ratio. The situation 
in the J /if} spectrum is somewhat worse. There, Davies et al. []23| find p = 1.2(2), while the 
experimental value is 0.48(1). Thus it appears that higher order corrections in the inverse 
quark mass may become substantial at the charm quark scale. Clearly further lattice work 
is required to settle this issue. The effects of unquenching may also be studied theoretically 
(although this will necessarily be model dependent). For example, Eichten et al. p|J have 



used the Wigner-Weisskopf method in a vector confinement model to study the importance 



of virtual meson decays. Geiger and Isgur |25| have also studied such effects on the structure 
of the nucleon using the 3 Po model to incorporate meson loops. 

The methodology which we have adopted here - using the heavy quark expansion of the 
Coulomb gauge QCD Hamiltonian to identify pertinent matrix elements and then using the 
(extended) Flux Tube Model to calculate these - is potentially very useful. For example, 
we may calculate the spin independent shifts to the static potential in precisely the same 
manner. This should prove interesting since we expect that they may not be obtained 
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through the nonrelativistic reduction of an assumed kernel. Rather we will see for the first 
time evidence of the dynamical nature of confinement in the effective 1/m 2 structure. Spin- 



orbit splittings in baryons have been something of a mystery for a long time P6fl . These 
should be accessible with the same techniques used here - assuming that one can find a 
satisfactory quantum flux tube model of gluonic excitations in baryons. It should also be 
possible to examine the strong decays of heavy quarkonia with this approach. This is of 
interest since these processes are poorly understood and are of central importance to many 
phenomena of current interest. Finally, most current models of hadrons do not contain 
dynamical (nonperturbative) gluonic degrees of freedom. It should, however, be possible 
to include them in a way which is consistent with the Flux Tube Model. We are currently 



investigating this possibility in the context of the Dynamical Quark Model [15 
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